arXiv: 1503.05935v2 [astro-ph.HE] 11 Jan 2016 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 12 January 2016 


(MN style file v2.2) 


Inverse-Compton drag on a Highly Magnetized GRB jet in 
Stellar Envelope 

Chiara Ceccobello^’^, Pawan Kumar^ * 

^ “Anton Pannekoek” Instituut voor Sterrekunde, Universiteit van Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands 
^ Dipartimento di Fisica, Universita di Ferrara, via Saragat 1, 441^^, Ferrara, Italy 
^Department of Astronomy, University of Texas at Austin, TX 78712, USA 


Accepted 2015 March 2. Received 2014 July 16 


ABSTRACT 

The collimation and evolution of relativistic outflows in 7-ray bursts (GRBs) are 
determined by their interaction with the stellar envelope through which they travel 
before reaching the much larger distance where the energy is dissipated and 7-rays 
are produced. We consider the case of a Poynting flux dominated relativistic outflow 
and show that it suffers strong inverse-Compton (IC) scattering drag near the stellar 
surface and the jet is slowed down to sub-relativistic speed if its initial magnetization 
parameter (cto) is larger than about 10®. If the temperature of the cocoon surrounding 
the jet were to be larger than about 10 keV, then an optically thick layer of electrons 
and positrons forms at the interface of the cocoon and the jet, and one might expect 
this pair screen to protect the interior of the jet from IC drag. However, the pair screen 
turns out to be ephemeral, and instead of shielding the jet it speeds up the IC drag on 
it. Although a high uq jet might not survive its passage through the star, a fraction 
of its energy is converted to 1-100 MeV radiation that escapes the star and appears 
as a bright flash lasting for about 10 s. 

Key words: (stars:) gamma-ray bursts: general - stars: jet - stars: magnetic field - 
scattering 


1 INTRODUCTION 


Long Gamma Ray Bursts (long GRBs) are explosions resulting from the core collapse of massive stars at the e nd of their nuclea r 
burning life cycle . The amount of energy produc e d in these exp l osion s is estimated to be 1 0‘^^-10®^ erg, e.g. (199£ h 


12001 




Cenko et al 


Dai^et_^]j[2001 Panaltescu fc Kumail ( 2001 1: Berger et al. (2003); Gurran et al. 1 2008h : Liang et al. 1 2008h : 


_ _ _ , _ _. . , _ . , _ ... . , R.acusin et al 

i 20inh . The collapse of the core of a GRB progenitor produces either a blac k hole or a neutron star 


and in either case t he cen t ral compact object is belieyed to be rapidly rotating (for a reyiew Piran ( lOOfll ): Meszaros ( 200fil ) : 
Woosley fc Bloom i 20061 ) : Gehrels et al. ( 20091) ). As in other astrophysical sources such as active galactic nuclei, micro¬ 
quasars, pulsars and SGRs (soft-gamma-ray repeaters) a rapidly rotating black hole, or a magnetar, is expected to pr o duce a 
relativistic bipolar je t which then int e ract with the ambient medium (e.g. BucciantinLet_aL ( 20081 . 20091 ) Falcke et al. (2004) 


Markoff et al.1 l|20oii ) ; Markofil ( 2O10l ) , Naravan fc McGlintockI ( 20081 ) , Yuan fc^a^w^ l|20ld )~ 


In the scenario where the central engine of a GRB is a rapidly rotating magnetar, 
genera t ed by the strong ma gnetic field with an initial jet magnetization parametei0, fro, 


a Poynting-flux dominated jet is 
of order ~ 10® jThompson et al 


( 2004h : iMetzger et all (120071) ). The magnetization parameter increases as the neutrino driven baryonic mass loss r ate at the 
surface of the neutron star decreases on de-leptonization timescale of about half a minute ( Metzger et all . I2OI1I ). In fact. 


the increase to the m agnetization parameter can be rather dramatic with cro ~ 10® as the neutrino luminosity winds down 
( Metzger et al . 201lh . These authors associate the transition to high cro with the end of the prompt 7 -r ay phase and the 


steep d ecline of X-ray after gl ow that is seen f or a l arge fraction of bursts detected by the Swift s atellite ( Taglia|erri__et_al 


( 2 OO 5 I) : O’Brien et al. ( 20061 ): Willingale et al. ( 2010| )). The reason for this association, according to Metzger et al. ~ (l2011 ). i 


* E-mail: c.ceccobello@uva.nl; pk@astro.as.utexas.edu 

® Magnetization parameter is defined as the ratio of Poynting flux and particle kinetic energy flux carried by the jet. 
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2 Ceccobello and Kumar 


that the acceleration and dissipation are very inefficient processes for high ao jets. It shonld be noted that highly magnetized 
jets are not limited to the magnetar model, but could also be produced when the GRB central engine is an accreti ng black 
hole as the mechanism for launching of the jet might be the Blandford-Znajek process (IBlandford fc Znaiekl . 119771 ). In this 
paper we address the question regarding the survival of a highly magnetized jet (ao > 10"*) as it propagates through the GRB 
progenitor star and is exposed an intense radiation field that ca n penetrate all t he way to jet ax i s and c an drag it down. IC 
drag has been consider e d by a number of people for AG N jets e.g. Phinnev I IQStI I: Melia fc Konig] I IQSflh : Sikora et al. I IQQfil h 
Ghisellini et al. ( 2005ll : Ghisellini fc Tavecchio ( 20101 1. However, the GRB jets are different in that they are highly opaque 


throughout much of the region where they undergo acceleration (as opposed to AGN jets which are transparent throughout 
their entire length), and thus the IC drag on GRB jets due to an external field can be treated independently of the acceleration 
mechanism. 

A jet propagating through the GRB progenitor star envelope creates a hot cocoon within which the jet is enclosed by the time 
it reaches the stellar surface. We describe the interaction of such highly magnetized jets with photons from the hot cocoon 
and how that affects the jet propagation through the envelope of the progenitor star. In particular, we analyze the role of the 
inverse Compton drag in braking the relativistic outflow. 

The organization of this paper is as follows. In the next section we calculate how far inside the jet can X-ray photons from 
the cocoon penetrate the jet, and at the distance from the center when the jet becomes transparent. §3 provides an estimate 
for IC drag on a high magnetization parameter jet due to scattering of X-ray photons from the cocoon surrounding it, but 
ignoring the creation of pairs and possible shielding of the jet provided by these particles. In §4 we consider whether the 
core of the jet could be shielded by electron-positron pairs produced either at the base of the jet or at the interface of the jet 
and cocoon by collision of X-ray photons with 7 -ray photons that are arise when X-ray photons are inverse-Compton (IC) 
scattered by electrons in the jet. 


2 PHOTOSPHERIC RADIUS FOR POYNTING FLUX DOMINATED JETS 


Let us consider a relativistic Poynting-dominated jet composed by a mixture of baryons, leptons and photons. As seen in the 
frame comoving with the flow, baryons and leptons are thermally distributed and have the same number density. 

The total isotropic equivalent luminosity of a Poynting jet, which has significant thermal energy, and where mass flux is 
dominated by baryons, can be written as 


L = -KR^e^R) 


f 2^.2 ,4 ,^2 , 

lUpUpVa 1 -I- V - - - 


( 1 ) 


where T and v are jet Lorentz factor and speed, n'p and B' are jet comoving frame proton density and magnetic held strength 
Uy is the energy density in photons in jet comoving frame, and 0j{R) is the jet half-opening angle when it is at radius R. 
Using the mass conservation equation and the dehnition of the magnetization parameter, which are the following 

M± = -KO^R^mpnpTv 

B’2 


AirmpUpC^ ’ 


( 2 ) 

(3) 


we obtain 


where 


L = M±c=r (1 + e + cr), 




4uL 


Smpn'pC^ 


(4) 


(5) 


is the ratio of thermal energy and baryon rest-mass energy densities. From equations ©-O we calculate proton number 
density in the comoving frame 

L 1 L 


n'p{R) 


( 6 ) 


ro(l -I- ^0 + cro)c^ 'KOjR^mpcT TrO^R^mpC^iao -I- ^o)r ’ 
where Fq, ^0 and ao are values at the base of the jet at radius Rq, and the last part of the above equation is obtained by 
assuming that Fq ~ 1 and uo -I- ^0 S> 1. 

While the jet is inside the star it is collimated by the pressure of the cocoon and the ambient stellar medium, and we 
take its Lorentz factor to increase with radius as 


r 

Ro 


(7) 


as long as F < (^0 + uo). The index a in the above equation can be shown to be re lated to the pressure (p) stratification of 
GRB progenitor star; for p oc r““, a = a/4 as long as a < 2 ( Kumar fc Zhan3 . 2014 1. The pressure in the He-envelope of the 
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IC drag on Highly Magnetized GRB jets 3 


GRB progenitor star declines as ~ r an d therefore we expect a ~ 0.5 for a Poynting jet. The transverse size of the jet 
increases with radius in the same way as F (IKumar fc Zhang . bPlJ l. i.e. 


~ Ro 


r 

Ro 


and therefore the jet angle is given by 


e(r 


Ro 

r 


( 8 ) 


(9) 


The photospheric radius where the Thompson optical depth of the jet in the jet-longitudinal direction is unity, is deter¬ 
mined from 

4a-l 


TT 


/ dr / , 

— aTn^{r)T ■. 


(jj'L 


Ro 

”r 


( 10 ) 


27r(4a — l)Romp(?{(Jo + Co) 

where we made use of equation for electron number density (which is assumed to be equal to proton number density), and 
equations 0 & ® to substitute for F and 6j respectively. Thus the photospheric radius is 

1 

(7tL 


Rph — Rq 


(2.4 X 10®cm) L5oao,6 (^20j,-i) 


( 11 ) 


27r(4a — l)mpC^{ao -|- Co)Ro_ 

So in case of initially large magnetization parameter (> 10"*), photons can escape from the jet well before the jet reaches 
the stellar surface. We note that the maximum Lorentz factor of a thermal fireball with no = 0 is obtained when the jet 
acceleration continues out to the photospheric radius Rph and is given by Vmax ~ [aTL/{2n{4:a — as long 

as Rph < R*; for a = 1, Fmax ~ 5.4 x 10^Rj^^. 

Photons from the cocoon surrounding the jet cannot penetrate very far inside the jet at radius Rph because of the much 
larger optical depth in the transverse direction. To estimate the inverse-Compton (IC) drag of the jet due to scattering of 
X-ray photons from the cocoon by electrons in the jet we first determine the radius where the jet becomes transparent in the 
transverse direction. 


2.1 Jet transparency radins in the transverse direction 

The optical depth of the jet for photons of frequency Vc traveling in a direction perpendicular to the jet axis, including 
Klein-Nishina correction to the scattering cross-section, is 


T±{r) 


aTnp{r)r9jr 
1 + hucF / (meC^) 


1 -h 2 In 1 -f 


2hi>cX' 

rrieC^ 


(jtL 


nOjiripC^lao -I- Co)’'[l + hvcT /(meC^)] ’ 


( 12 ) 


where we used equation m for Up. The radius where the jet becomes transparent to photons moving in the transverse direction 


Rph,± « (5 X lO^^cm) 


1^500-0,6 


(13) 


A more general situation is where photons are traveling at an angle 6-y with respect to the jet axis. We calculate the 
optical depth for these photons to travel from the interface of jet and the cocoon to the jet axis. 

Consider a photon traveling at an angle 6~f with respect to jet axis and electrons moving in the radial direction (see 
Fig.0. The optical depth for a photon of frequency Vc to scatter off electrons along its trajectory starting from cocoon-jet 
interface to the jet axis is 


T{e-, 


F" aTnp{r)[l-vcos{e + e^)/c] 

/o [l + hvGmeC^\s\Ti{9+ 0^) 


(14) 


where rip = Trip is electron density in star rest-frame, 

u'R9 + 9j) = UcT [1 — vcos{9 9^)/c\ , (15) 

is photon frequency in electron rest frame, and we have assumed that there is one electron per proton in the jet; Klein-Nishina 
correction to Thompson cross-section is included in eq. m- The photon trajectory in the polar coordinate is described by 
(see Fig. 0 

r sin(0-I-0.^) = 7?sin(0j-I-0.y). (16) 

Using 


,(r) = 


TT9jr^mpC^{ao -f ^o) 


= (7 X 10^'^cm"®) 


7/50 


^?,-l»'lo('^0.6 +C 0 . 6 )’ 


we find 


(JtL 


Tr9jmpC^{ao + fo) Jo 


d9 


[1 — vcos(9 + 9~,)/c\ 


r[l -I- hvR9 + 9^)/meC^\ sin(0 -|- 9ry 


(17) 


(18) 
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4 Ceccobello and Kumar 



Figure 1. Schematic sketch of a photon trajectory (red line) from the hot cocoon passing through the jet axis (dotted line). The angle 
between the photon initial direction of motion and the jet axis is labeled as 6^ and the jet opening angle is 9j. R is the radial distance 
between the center of explosion and the injection point where the cocoon photon enters the jet and r is the radial coordinate. Although 
the jet is shown to be conical in this representation, its transverse radius increases more slowly than r while it is confined by the pressure 
of the star/cocoon. 


With use of equation m for r this reduces to 


21 -1 


T(r) 


arL [l + 

■K6^mpC^{ao + ^o)rsm{9j + 6j) 


\9j — V {sin(0j + 9j) — sinS^} /c] . 


(19) 


The GRB jet Lorentz factor near the surface of the progenitor star is expected to be much larger than 9- ^. In that case 


for 9^ 9j equation mi) reduces to 


T Rri 1.5L5o(cro,6 + ^0,6) [l + hu'c{6^)/niec] 


21 -1 


( 20 ) 


and for 9-y ';$> 9j the expression for optical depth reduces to equation m- The jet optical depth in the transverse direction at 
a given radius for these two limiting cases differs by a factor ~ 10. 


3 INVERSE COMPTON DRAG FOR POYNTING FLUX DOMINATED JETS 

We calculate the inverse-Compton loss for an electron when it is exposed to a beam of photons moving at an angle 9^ 
with respect to electron velocity. The electron moves with the jet and therefore its Lorentz factor and velocity are T and v 
respectively. Let us consider the specific intensity of the photon beam from the cocoon in the rest frame of the star to be 
Iv{9ry). The transformations of photon frequency, specific intensity, and angle from star-rest frame to jet comoving frame are 
given by 

v'— uT(l — V cos9-y / c), lli{9') = /v)^, sin9'= {u/v')sm9-./. (21) 

The rate of loss of energy for an electron due to inverse-Compton scatterings is proportional to photon energy flux in its 
comoving frame. Using the above transformations the comoving flux can be shown to be proportional to 0^r^(l —u cos 9jIc)‘^ oc 
for rS-, > 1. Thus, the IC drag on electrons increases very rapidly with increasing 9j, and that nearly compensates for 
the smaller flux at the jet axis due to larger optical depth for larger 9-y. So from here on we specialize to 9.y ~ 1. Consider a 
photon of frequency Vc scattered by an electron of Lorentz factor T. In the comoving frame, the scattered photon energy is 

hv'r 


hu'i^ = 


1 -I- (1 — COS9y) 

mec-^ ^ ' 

The average energy of the scattered photon in the Klein-Nishina limit can be calculated as 

I ^^7 dn " 


(hvic) = 


Jdn'^ 


( 22 ) 


(23) 


where the differential Klein-Nishina cross section is taken from jBlumenthal fc Gould . 197Clh . When 2 < hvc^j < T 


the average energy of the scattered photon can be shown to be 


{hvic) 


2mecV 


1-1-2 ln[l -I- 2hUcT / (rrieC^)] 


SrrieC 


2hPrT 


+ 1 - 


hUrV 


In 1-h 


hi/cT 

2m.c^ 


(24) 
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IC drag on Highly Magnetized GRB jets 5 


The equation for IC drag is 


dmpCTo 


Fc{t) 


U kn hVic 5 


(25) 


dt huc 

where is number of electrons and positrons per proton, Fc is cocoon thermal flux given by equation (IA.23II . and the 
electron-photon scattering cross-section as a function of dimensionless photon energy in electron rest frame (hVcF is 


Ok-niJlVcF 


Sctt 


1 -f 2 In 1 -f 


2hVcT 


ra.c‘‘ 


(26) 


l&hVcF / (nisC^) 

In deriving equation ( 1 ^ we assumed that electrons, protons and electro-magnetic fields are coupled and move together; 
charge particles are coupled to the jet EM field via E x B drift force which ensures that particles move with EM fields, and 
thus the total energy per proton in the jet is ~ ampC^T. 

The IC cooling time in the stellar rest frame for an electron in the jet, at radius r and time after the formation of 
cocoon, follows from equations (1251) and (1261) : 


^kn SmpC^ao{huc/meC^)‘^ 

~ SarUFcitr) 


3meC^ 

2hUcT 


+ 1 - 


rUeC 

hVnT 


In 1 -f 


hUcT 

2meC^ 


(27) 


Substituting for hvc = SksTc lea. IA.161) and Fc leg. IA.2^ we obtain the IC cooling time in Klein-Nishina regime to be 


die ~ (0.8s)exp(Tx)ty^o-o,6?7c,i''^C±^ 


(28) 


where rje is the terminal Lorentz factor of the cocoon. The IC cooling time when the scattering is in the Thompson regime, 
i.e. hVcT rrieC?, is given by 


(5x10 ^s) exp(Tx)ty^r 2 ^(Tq, 


6»7, 




(29) 


The dynamical time at the stellar surface is i?*/c ~ 3 s. The IC cooling time is smaller than the dynamical time as long as 
fro ^ 10®. Equation ([271) can be solved to find the radius, Rcc, where the IC drag timescale at the jet-axis {tic) is comparable 
to the dynamical timescale of the jet {tdyn ~ r/c), namely 

8mpRaotR^{huc/irieR)'^ ^ Rcc 


8aT(TBT^tF^ C± 


(30) 


where tr is the time when IC cooling is considered — it is in general larger than the dynamical time since cocoon formation 
begins with the launch of the relativistic jet and jet duration should exceed R^/c in order for the jet to break through the 
stellar surface — and tfs (given by eq. IA.22I) in the mean time in between scatterings of a photon while inside the cocoon; 
the equation is valid for hv^/{mc(?) 1 . 

Making use of equation dm) for Tx we can rewrite the above equation for Rcc in a more explicit form: 

3(TTCr Cx Rcc 


exp 


(jtL 


n6jmpC^{ao + ^o)I?cc[l + hUcT/{meC^)] 
Equation (1311) can be rewritten with the use of equation (|7|) for T 


SrUpC^aotR^ {hUc/meC^)^ 


^a/Rcc ^ KR^c 


where 


A = 


(JtL 


'K6jmpC^{ao -f ^o)[l -f hVcL/ {mec'^)] ’ 


K = 


3aTCJBTH)i\± 


(31) 


(32) 


(33) 


SnipC^aotF^ (hi'c/mcC^)'^ 

The equation for the Compton cooling radius Rcc is a transcendental equation which can be solved perturbatively after we 
transform it into the following logarithmic transcendental equation 

A 


Rcc 


= log)!!") -I- log(7?cc) 


At the zero-th order, log(7?cc) term on the right side of the above equation is neglected and, the solution is 

which, when substituted back into (IMl) . provides the first order solution for the Compton cooling radius 

A 


Rcc = 


(34) 


(35) 


(36) 


The estimate for Rcc using equation 


log(A') -I- [log(A) - log(log(AA:))] 
and the exact solution of eq. [32]for parameters Tso = l,6j = 0.1, Tc ,7 = l,r 2 = 


1, fro ,6 = 1, and C± — 1 nre Rcc = 7.5 x lO^^cm and Rcc = 7.0 x lO^^cm respectively. 
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6 Ceccobello and Kumar 


4 SHIELDING JETS FROM IC DRAG BY CREATION OF ELECTRONS AND POSITRONS 

Thus far we have assumed that the jet energy is carried outward by magnetic fields, protons and electron-positron pairs. 
However, we have not estimated the number of that might have been produced in the hot plasma at the base of the jet, or 
pairs that might be created in the collision of inverse-Compton scattered photons. The presence of these pairs could shield the 
inner part of the jet from IC drag due to thermal photons from the cocoon. In the next subsection we take up the calculation 
of thermal that owe their existence to the initial hot plasma at the base of the jet, and show that their number density 
is too small far away from the jet launching site to be able to shield the jet. In section 14.21 we provide an estimate of the 
density of pairs generated when thermal photons from the cocoon collide with photons that are IC scattered by in the jet. 
This process is shown to be effective in shielding the jet for a while but eventually pairs annihilate and pair screen disappears 
exposing the jet-core to severe IC drag f t|4.2l) . 


4.1 Thermal pairs and shielding of Poynting jets 


The number density of thermal pairs at any r is given by the standard thermal distribution formula corresponding to the local 
temperature of the jet as long as the annihilation time is less than the dynamical tim^. The radius where the two time 
scales become equal, Rfreeze, is the freeze-out radius for pairs. Beyond this radius the total number of does not change 
barring the dissipation of jet kinetic, or magnetic, energy and using that to create new pairs; non-thermal pair creation will 
be taken up in ^4.21 We calculate the number density of pairs at Rfreeze and show that to be much smaller than the density 
of protons. Therefore, thermal pairs are unimportant for shielding the jet. 

The temperature at the base of a Poynting jet is 

1 1/4 


fcsTo ~ ks 




47ri?ocrs(o-o -|- ^o) 


(41keV) Ro 


l/4p-l/2 


(37) 


which is considerably smaller than the temperature for a thermal fireball with ao < 1; L*®" = ALjO^ is isotropic equivalent of 
jet luminosity {6j ~ 1 at the jet-base). 

Considering the conservation of entropy in a shell of plasma as it moves to larger radius with the jet we find the decrease 
of comoving frame temperature with r 


T'{r) ~To[ — 


2q!/3 


n-1/3 




(38) 


where we have made use of equation © for the transverse size of the jet and equation © for r; as noted above equation 
Q ~ 0.5 in the helium-envelope of CRB progenitor star. 

The cross-section for pair annihilation when the average thermal speed of is u± is (tt/(u±/c). The annihilation time 
for a positron in jet comoving frame, given the density of electrons to be n±/2, is therefore 

^ (39) 


f 

^ann 

CFtW^C 

The pair annihilation ceases, and their total number freezes, at a radius where t), 
given by 

2 r 


R 


■freeze ' 


and the pair density at Rfreeze is 


n'±(Rf 

r 


aTn'f- ’ 


2V{Rfr 


r/cV{r). Thus the freeze-out radius is 

(40) 


(TT^fr 


To determine the pair freeze-out radius we substitute for thermal pair density, i.e. the following equation 

2{2-KkBmeT'f/'^ _^ _ _2 . 


n± 




■ exp {—mec/ksT'^ 


and r dependence of P and T' into (|40l): 


exp 


5.9 X 10^ (ILY 
To \Ro) 


p t-t3/2 
ttpiQ 

6.2 X 10® 


5ck/2-1 


Let’s define 


C = 


5.9 X 10® 

To ’ 


and 


3/2 


6.2 X 108’ 


(41) 

(42) 

(43) 

(44) 


^ A more precise statement is that number density of pairs at any given radius is given by the balance between creation and annihilation 
rates. However, it can be shown that the assumption of thermal equilibrium is approximately valid as long as annihilation time is short 
compared with the expansion time. 
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IC drag on Highly Magnetized GRB jets 7 


and rewrite equation (US as 


R 


c ( ;^ ) = log(-D) + 1 o« - 1) log ( — 


Ro 


R 


(45) 


which is easier to solve analytically. Neglecting, at first, the logi?-term on the right hand side, we find R = ((logD)/C)^^“. 
Substituting that back into (1451) . the approximate analytical solution for the freeze-out radius is found to be 


f tifr 


V Ro 




\og{D) + 


(log(log(D)) - log(A)) 


If O. 


(46) 


The results for the freeze out radius obtained as exact numerical solution of equation (1431) are in good agreement with the 
above approximated expression. 

The freeze-out radius has a weak dependence on a, and for ag = 10® and Tg = 41 keV, pair density freezes out fairly 
close to the jet launching site. 

The temperature at freeze-out can be calculated using equation 1311), and it can be shown to be ~ 10 that is 

almost independent of various parameters. The isotropic equivalent luminosity carried by pairs for r > Rfreeze is 


; 47rJ?f 


3 / 7^2 
grrieC nj_i 


‘^'R^Rfreeze'^eC T jc'J' 


47ri?omeC^ / Rf^ 


where we made use of equation m for pair density at Rfreeze- or 


L ISO 

± 

Rise 


5 X 10 


-16 Ro,7 


L ISO 

52 


reeze 

-ee. 

O T ^SO-^/ ^ p 

H, 


Rg 


l+Sct 


1/4 1/2 -1/4 


(3c.-l-l)/a 


(47) 


(48) 


e)Ll7 ^Rg,r 


(49) 


The jet kinetic luminosity at r is Ljair), and we see from the above equation that thermal pairs are insignificant carriers of 
jet kinetic luminosity — most of the kinetic luminosity is being carried by protons. The ratio of pair to proton number 
density above the freeze-out radius is given by 

n'p 

where cr(Rf reeze) ~ o-g/T{Rfreeze) IS the magnetization parameter at r = Rfr 
back into equation (|49l) it becomes 

nj 6 f Rq 

~ 0 CTO,6 ( ^ 

'^p \ ^ ^jree: 

Equation (I50|) shows that thermal pairs are too small in number to affect the propagation of photons into the jet, and 
hence they cannot shield the jet from the severe IC drag. 


e. Substituting this expression for o-(Rfreeze) 


1^2 i?0,7 2L: 




(3a+l)/a 


0.7 J 


(50) 


4.2 Pair creation due to photon collisions and shielding of jet from IC drag 


The calculation in previous sections ignored the possibility that thermal photons IC scattered by the jet might have sufficient 
energy for electron-positron pair creation. These IC scattered photons could give rise to an optically thick layer of on the 
side wall of the jet that is in contact with the cocoon, and this could potentially shield the interior of the jet from IC drag. 
We investigate that possibility in this sub-section. 

Let us consider the mean frequency of thermal photons in the cocoon to be Ve, and after colliding with an electron in 
the jet with Lorentz factor T the frequency increases to Vie (these frequencies are in the rest frame of the star). The scattered 
photon travels within an angle T”*^ of the electron’s velocity vector due to relativistic beaming, and hence its chances of 
undergoing a second collision with another electron is drastically reduced since the photon is moving in nearly the same 
direction as electrons in the jet. 

The condition for pa ir production when averaged over the angle between colliding photons is: {hvic){hve) > e.g. 

Gould fc Schreden ( 19671 1. Since hvic ~ me(?Tl2 (eq. 12411 for large angle collisions in the Klein-Nishina regime - i.e. when 


hVcV > rUeC? - the pair production condition becomes hUeT ^imeR. Thus, for a given cocoon temperature (Tc) at a certain 
radius where jet Lorentz factor becomes larger than the following critical value 


4meC 

SksTc 


(51) 


pair production at the interface of the jet and cocoon begins. For cocoon temperature of ~ 10 keV, Terit ~ 70. The Lorentz 
factor of magnetic jets increases with radius as ~ and therefore T ~ 10^ can be attained near the stellar surface where 
the jet also becomes transparent in the transverse direction (in the absence of e^) for ag > 10®; the Lorentz factor of a 
thermally driven jet increases more rapidly with r, and it too is likely to develop an layer surrounding it before reaching 
the stellar surface. 

For the remainder of this section we assume that the condition for pair production is satisfied, and describe the effect 
that has on jet structure and dynamics. 
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= (cr±n^) ^ ~ (2 X _i) 


The distance traveled by a high energy IC photon — which is moving almost parallel to the jet axis (within an angle T ^ 
to be precise) — before it is turned into as a result of collision with a thermal photon from the cocoon is 

(52) 

where cr± is the cross-section for 7 - f 7 ^ e~ - | - e"*" ( the maximum value for a± is 2.5 x 10“^® cm^ at photon energy 1.4 times 
the threshold value given above e.g. ISvenssonI 1 1982l ~l. 

hVcC 

is the number density of thermal photons at the interface of the cocoon and the jet at time t (in seconds) in star rest-frame, 
t'c = SksTc and Fc{t) are given by equations (IA.16I) and (IA.23I) . We see from equation (1521) that high energy IC photons do not 
travel very far from their place of creation before undergoing pair production. The newly produced e®^s have thermal Lorentz 
factor less than ~2 in jet comoving fram^ but they cool rapidly via the synchrotron process on a timescale much smaller 
than the dynamical time; magnetic field in the jet comoving frame is B' = {AL/O^T^r'^cY^'^ = (1.2 x lO®G)LgQ^/[0j,_ir2rii], 
and hence the synchrotron cooling time in jet comoving frame is ~ (5 x 10~®s) [ 6 j^-ir 2 riiYL^Q . 

The rate of pair production per electron is approximately equal to the number of photons it scatters per unit time, i.e. 


n± 

n. 


(TT'fl-yC 


(3 X 


(54) 


1 -I- hvcT/im^cY 

The second equality is obtained for the case where thermal photons from the cocoon are scattered by electrons in the jet 
in Klein-Nishina regime, i.e. hVcB = 'iksTcTj 1, and in that case we find the surprising result that the 
rate of pair production per electron depends only on the Lorentz factors of the jet and cocoor0. Once pair production starts 
it proceeds rapidly since newly produced also IC scatter thermal photons which have sufficient energy for more pair 
production. The number of pairs produced per proton is expected to be of order mpT/{meVcrit)', the number is limited by 
pair production rate (which drops to zero when jet LF falls below Tcrit and IC scattered photons no longer have sufficient 
energy for pair production) and annihilation rate. Thus, pair production saturates on a timescale 

tsat ~ (3 X 10“®s)t^^^ r2 In [mpT/meFcrit] . (55) 

The physical thickness of the pair screen corresponding to optical depth unity, A±(t± = 1 ), is given by 


A±(t± = 1) = 


(me/rup) 


(1.2 X 10®cm) 


^»,ii(no,6 + ^O.60|,-l) 


(56) 


UTUp (T/Pcrit) L50 

where rip is proton density in the jet at r = /?* which is given by equation GZl), and ^0 is thermal energy per proton divided 
by rripC^ at the base of the jet. 

We note that as long as the Lorentz factor of is larger than Tcrit, 7 -rays from pair annihilation are converted back to 
within a short distance of a few meters (eq. 1521) due to collisions with thermal photons from the cocoon. Pair production 
process is terminated when the Lorentz factor of the screen falls below Tcrit and its optical depth becomes so large that X-ray 
photons from the cocoon are prevented from entering the faster moving interior of the jet where F > Tcrit- 

Let us consider that the optical depth of layer surrounding the jet at radius r is r±(r). The flux of photons from the 
cocoon that crosses this layer is 


U., it,T±) 


Fcjt) exp{-T±) 
iksTc 


(5.5 X 10®®cm¬ 


's )exp(-r±) 




(57) 


where we made use of equations (IA.I 6 II and (IA.23II . Pair production can proceed in the interior to the e^ screen as long as 
the distance traveled by a IC scattered photon before it collides with a thermal photon (eg. 1521) is smaller than i?*, i.e. 


Aj^( 2 x lO^cm) e"± (Lso/Sj.-i) 


1 / ^^*Vcp^<R* 


(58) 


Moreover, the time it takes for a photon from the cocoon to cross this layer should also be less than i7*/c. Thus, we hud that 
the optical depth of the screen is 

r± ~ 20 + log R:p{\Ko/e,,-iY^Wc!i] , (59) 

and the physical thickness of the pair screen is max{20A±, i7,/r}, where A± is given by equation (1561) : pair screen thickness 
is Rt/T when IC photons — moving at an angle with respect to the jet axis — travel a distance i?, before colliding with 
thermal photons from the cocoon. 

Electrons and positrons in the pair screen continue to scatter X-ray photons from the cocoon and the resulting drag slows 
down the jet below Tcrit (hg-Ill)- The time it takes for F to fall below Tcrit is of order 10~^s exp(rx)/C± lea. 1281) . When the 


® Newly produced pairs are swept up by the magnetic field in the jet and forced to move with the outflow. 

The time dependence for pair production rate of (eg. 1541) is due entirely to the cocoon thermal flux Fc- 
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Figure 2. Schematic view of the layered jet structure while it has been braked by the Inverse Compton interaction with the hot photons 
from the cocoon. If the magnetization parameter a is larger than 10^, the jet build up a self-shielding e'^/e“-pair screen created by 
photons from the cocoon and IC-scattered photons. 


Lorentz factor of pairs falls below Tcrit feo. I51[). 7 -rays produced in pair annihilation no longer have sufficient energy for pair 
creation, and at that time the pair screen begins to evaporate. 

The time (measured in the rest frame of the star) for a positron to run into an electron and annihilate is 


te±,an{r) = 






(60) 


_>^^n±v± (TTn±c ' ' C±Ll2° 

where = gt l(v±lc) is annihilation cross-section, n± is the average relative speed between electrons and positrons in 

jet comoving frame, and n± ~ C±np is the pair density; for ((± ~ {irip/2ms)(T/Tcrit), i-e. the maximum possible number of 
pairs per protons when F > Tcrit, the annihilation time is ~ 0.3 s. The ratio of annihilation and dynamical times at radius r 
is 


tc±,ar{r) ^ ^ ^ ^ ^Q2 ririi(ao 6 +^o,6) ^ 
tdyn 

which is smaller than 1 for r < 10 ®cm even for = 1 , and that means that any pairs produced at the base of the jet or at any 
radius smaller than 10®cm cannot survive as the jet moves to larger radii. Hence, only an ongoing process of pair formation 
can support (C± > 1 . 

The Lorentz factor of pair screen continues to decrease due to IC drag and that causes the annihilation time — which 
scales as F® (eq. lOOlfl — to decrease rapidly. Since the timescale for the creation of pair screen is smaller than the IC drag 
time which is much smaller than the pair annihilation time, soon after the pair screen forms its Lorentz factor decreases 
rapidly to order unity (in 1 /rs or less - eg. 1291) . pairs annihilate and screen evaporates on timescale of order 1 ms fea. 1601) . 
Once the optical depth of the screen decreases, photons from the cocoon pass through it, and the process of formation of e* 
and IC drag progresses deeper inside the jet. At any given radius this process is only terminated at a distance from jet axis 
where the Thompson optical depth out to the jet-cocoon boundary, for ~ 1, is of order 10. Hence, as we get closer to the 
stellar surface we find an increasingly larger fraction of the jet to have gone through the process of forming pair screen, e* 
annihilation, and IC drag. If the jet with = 1 becomes transparent in the transverse direction near the stellar surface - as 
it in fact does for go > 10® according to eq. [13]- then self-generated pair screen is too short lived to protect it from IC drag. 
It should be noted that the transient nature of the pair screen in fact speeds up the process of jet IC drag because the drag 
time is proportional to (see eas. 12812^ . 


® For a time-independent relativitic outflow, particle density in star rest frame does not vary as we follow the flow-steamlines even though 
the Lorentz factor might increase or decrease. This is due to the conservation of particle flux and the fact that the speed is a constant c 
for a relativistic system. The particle density along a flow-line, of course, varies as F”^ in the comoving frame of the outflow. 
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Although very high-cr jets (ao > 10®) are unlikely to survive their passage through the GRB progenitor star and cocoon, 
they do not disappear without leaving an observational signature. Annihilation of pairs in the screen when the jet Lorentz 
factor decreases from Fcrit to of order unity produces photons of energy between 0.5 MeV and ~ meC^Fcrit ~ 50 MeV which 
can escape the jet in the longitudinal direction to arrive at the observer. The jet is slowed down by pair creation and IC 
drag, and these processes are about equally important for decelerating the jet. Hence, the observer-frame luminosity carried 
by pairs is of order the jet luminosity (Lj). The total energy carried by the photons resulting from pair annihilation is of the 
order of the energy of the pairs, and therefore, the total luminosity of pair annihilation photons, Lj, is of the order Lj. The 
duration of the annihilation pulse we expect to be of order the activity time of the central engine, which for long GRBs is 
typically around 5-100 s. 

The picture that emerges is that the outer layers of GRB jets are slowed down due to IC drag. However, inner regions 
continue to accelerate with r, and at some radius their Lorentz factor exceeds Fcrit- At that point a very rapid generation 
of ensues provided that the optical depth of the slower moving layer outside of this region is less than about 20. The 
IC drag slows down the pair screen rather quickly and then annihilate on a relatively short time scale of ~1 ms, and 
the formation of a new pair screen moves closer to the jet axis. This process continues until the jet becomes opaque in the 
transverse direction due to just the electrons associated with protons, i.e. for = 1. Jets of initial magnetization (cto) smaller 
than 10® are sufficiently opaque in transverse direction even when they rise above cocoon surface that they are essentially 
protected from IC drag. However, photons from the cocoon can penetrate a jet all the way to its axis when cto 10®, and the 
strong IC drag then slows down the outflow to sub-relativistic speed. It turns out that high-cr jets cannot escape this fate in 
spite of the e^ pair screen they create because these screens are rather short lived. 


5 CONCLUSIONS 

Relativistic jets in GRBs are surrounded by a hot cocoon of plasma that was created during the initial passage of the jet 
through the star when it shock heated the gas along its path and pushed it sideways to clear a cavity through the polar 
region of the GRB progenitor star. Thermal photons from this cocoon are scattered by electrons in the jet and that provides 
a strong drag force on the jet. Jets of initial magnetization parameter (ctq) smaller than about 10® are highly opaque in the 
transverse direction while traveling inside the GRB progenitor star, and thus they have a core region that is protected from 
this IC drag. The outer layers of this jet (about 20 Thompson optical depth thick), however, suffer IC drag and are slowed 
down considerably. 

Jets with CTO > 10® are transparent to photons from the cocoon near the stellar surface, and they are slowed down to 
sub-relativistic speeds due to IC drag. This is in spite of the fact that an optically thick layer of electrons and positrons forms 
at the interface of the cocoon and jet and that tries to protect the core of the jet from IC drag (fig.E}; these pairs are formed 
by the collisions of thermal photons from the cocoon with high energy photons that are produced when cocoon photons are 
IC scattered by the jet. However, the problem is that the pair screen itself slows down rapidly due to IC drag, and that causes 
pairs to annihilate and the e^-shield to evaporate rather quickly. Pair production then moves closer toward the jet axis, and 
the story is repeated until the entire jet is slowed down by the IC drag. 

The process of pair-screen formation and annihilation associated with a high-cr jet has observational consequences. For 
jets with Ctq > 10® we should see a pulse of high energy photons of energy between ~ 1 MeV and ~ nicC^Fcrit ~ 50 MeV with 
luminosity of order the Poynting jet luminosity. The duration of this pulse should be of order the central engine activity time. 

Even Poynting jets of ctq < 10® - which are highly opaque in the transverse direction — might suffer effects of IC drag 
indirectly. As the outer layers of these jets are slowed down by IC scatterings the resulting shear instabilities might slow down 
the inner regions as well. Moreover, if magnetic field lines thread the outer and the inner regions of the jet then slowing down 
of the outer part of the jet would get communicated to the inner region, and that could effect the entire jet. The details of 
this would depend on the magnetic field configuration and that is something that needs to be looked into. 

If a highly magnetized jet manages to escape IC drag while traveling inside the GRB progenitor star - for instance if it 
is encapsulated inside a highly opaque baryonic outflow to shield it from X-rays from the cocoon - it would be subjected to 
rapid dissipation due to charge starvation before reaching the deceleration radius (appendix B). 
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APPENDIX A: RADIATION FROM COCOON SURROUNDING RELATIVISTIC JET 


We describe in this appen dix a simplified, ana lytical, treatment of cocoon dynamics and radiation. This follows the work of 
Ramirez-Ruiz et al. ( 20021 1 and Matznei ( 20031 ) except for one thing and that is that we do not assum e that t he jet is conical 
in sha p e while inside the sta r. Numerous investigati o ns of relativistic jets, e . g. Lazzati &: BeeelmarJ (120051 1. iMorsonv et al. 
( 2007 1. Mizuta fc Alov ( 2009h . Mizuta fc loka ( 2013l l. Bromberg et al. J 2 OII . 2014h have shown that the cocoon created by 
the jet can be very effective in collimating it, and hence we take the jet opening angle to be a function of distance from the 
center. 

The energy in the cocoon (Ec) is of order the energy carried by the jet while it makes its way through the polar region 
of the GRB progenitor star, i.e. 

Ec ~ LRt,/vh, (A.l) 

where Vh is the average speed at which the jet head moves through the star. The jet head speed can be calculated from the 
conservation of momentum flux in the radial direction of the unshocked stellar gas as viewed from the rest frame of the jet 
head: 


P,c2(rj/4r|)« p^rlvl 


(A.2) 


where pj and pa are densities of the unshocked jet and the stellar envelope respectively, and Fj and Fjt are the Lorentz factors 
of the unshocked jet and the jet-head with respect to the unshocked star (the Lorentz factor of the unshocked jet with respect 
to jet-head is Tj/2rh)- Considering that the jet luminosity at the stellar surface can be written as 


and the mass of the swept up gas by the jet is 


we obtain 


L = -Ke]RlpH]c^, 

(A.3) 

rUc ~ ttO^PcRI, 

(A.4) 

2 \R,Ly^^ 

}^Vh ~ , 

crric 

(A.5) 


where Oj is jet opening angle which is a function of distance from the center due to collimation provided by the cocoon. We 
can simplify this expression further by substituting for L using equation dm 

4EhVh ~ cpc (A.6) 

where 


Vc = 


(A.7) 


is the terminal Lorentz factor of the cocoon plasma (provided that rjc > 1) after it escapes through the stellar surface and its 
thermal energy is converted to bulk kinetic energy. 

Therefore, the jet head speed is sub-relativistic when pc < 4 and is given by 

Vh ~ cpc/4. (A.8) 

For rjc > 4, the jet head speed is relativistic and its Lorentz factor is given by 

~ (Pc/4)'/^ (A.9) 

The expansion speed of the cocoon in the direction per pendicular to it s surface, Vc, is determined by equating the ram pressure 
with the thermal pressure inside the cocoon (pc), e.g. Matznerl ( 20031 ). 

(A.IO) 
(ATI) 


Vc = [PcIPaY^^- 


The average thermal pressure inside the cocoon is approximately 

E, 
P. 


3V 


where 

Vc ^ nRvhvl/3 ^ Rl{vc/vhf (A.12) 

is the volume of the cocoon at the time it emerges at the stellar surface. Combining equations (lA.lOl) and (lA.llll we find 

4 EcVh _ a2„..2 2, 


3p, 


HA A A /o 

^ ~ QiPcVhC /3. 


Substituting for Vh from equations (TOl) and dm we find 

Vc I for Pc < 4 


{Q]r]cl3f^^ for4<pc<6l-' 


(A.13) 


(A. 14) 
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12 Ceccobello and Kumar 


The thermal pressure of the cocoon can be obtained using equations (lA.lOII and (IA.14II and is given by 

L 


and its temperature is 


ejRic{2,n.y/^' 


keK = keiSpcKay^'^ ~ (24keV) 


-1/4 


fll/4 pl/2 1/8’ 


(A.15) 


(A.16) 


where aa is the radiation constant. We note that the cocoon temperature has a weak dependence on jet luminosity and angular 
size, and so it is unlikely to be larger than ~ 30 keV. 

The number density of thermal pairs at temperature Tc is given by 

2(27rfeBmeT,)3/2 


n± 


hi 


■ exp [—niec/kBTc) . 


(A.17) 


Therefore, for ksT^ = 20 keV, n± = 1.1 x lO^^cm and for 30 keV cocoon temperature n± = 1.1 x lO^^cm We next 
calculate the number of electrons associated with protons in the cocoon, and show that these exceed as long as ksT < 30 
keV. 

The average number density of electrons associated with baryons in the cocoon is 

Trie 3po t- ,, i„21_—3\ 7/50 


nipC^ric 


(1.5 X 10^ cm"'^) 


Or-iRlriVK 


3/2 ■ 


(A.18) 


The electron number density near the stellar surface, however, is smaller than the average density given above. The density 
at the stellar photosphere is ~ 1/{otHp), where 

Hp = Cl/g ~ ( 10 ®cm)T*, 57 ?J,iiM-i\ (A.19) 

is the density scale height, Cs is sound speed, T* is photospheric temperature, and is stellar mass in units of 10 Mq. 
Thus, the electron density at the photosphere is of order 1.5 x 10^® cm~®, and it increases with depth as ( 2 / 2 *)^^^'’'“'^^; where 
7 ~ 1.5 is the effective polytropic index that describes stratification near the stellar surface, and 

2 * =/f/p/( 7 -l). (A.20) 


Therefore, pair density in cocoon is larger than proton density at the photosphere as long as ksT > 15 keV, but at a depth 
of more than a few scale height below the photosphere the proton density exceeds n±. So the electron density in the cocoon 
as a function of radius can be written as 


n± + min 


arHf 


1 + 


{K - r) 


l/(7-l) 


, (1.5 X lO^^cm 




(A. 21 ) 


The time in between scattering for a photon (photon mean free time) in the cocoon is given by 


ifs 


(A.22) 


or c rie.c 

The thermal flux at the interface of the cocoon and the jet is dictated by diffusion of photons in cocoon, and is given by: 

Feit)^aBT*[tfs/{t + tfs)]^^^. (A.23) 


This expression is valid as long as t is less than the cocoon expansion time ~ Rt/vh- The flux at an optical depth r inside 
the jet is /c(t) exp(—r). 
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APPENDIX B: CHARGE STARVATION OF A POYNTING JET 


Let us consider a Poynting jet of isotropic equivalent luminosity L“‘° and magnetization parameter at its base of uq- The 
magnetic field in the jet is assumed to change direction on a length scale of Ib (in star rest frame) which corresponds to 
I'b = tsT in the jet comoving frame. The current required for supporting this non-zero curl is 

j' ~ B'cIGttI'b), (B.l) 

where 


B' = 


r- 


1/2 


(B.2) 


is magnetic field in jet comoving frame. 

The electron density in jet comoving frame for a Poynting jet of high magnetization parameter is obtained using equation 
© and is given by 


n'eir) 


(B.3) 


Anr^mpC^aoT ’ 

where is the number of per proton which we know from the discussion in section U should be of order unity for r ^ R^. 

The current required to support the jet magnetic field must be smaller than the maximum current that can be carried by 
charged particles in the jet, i.e. j' < qn'f,C- It follows from this requirement that beyond a certain radius, Res, the jet becomes 
charge starved, i.e. it does not have sufficient number of electrons to carry the required current. Using the above equations 
we find this radius to be 


Rc 


mpaoc° 


(1.5 X 10^'^ cm) 


riso^B,7T2 


(70,6 


(B.4) 


This should be compared with the deceleration radius, Rd, where the energy of the medium swept-up and shock heated by 
the jet is approximately half the total energy of the explosion: 


Rd 




noP^ J 


(1.2 X lO^^'cm) 


T-t^SO 

-^53 

noPi 


1/3 


(B.5) 


where is the isotropic equivalent of total energy carried by the jet, and no is the mean number density of protons in the 
circum-stellar medium of the GRB. Thus, a high magnetization jet that stops accelerating when it attains a Lorentz factor of 
~ will become charge starved near the deceleration radius and its magnetic field will dissipate rapidly. 

We note that if the jet were to start spreading in the radial direction at r < Res then it would never become charge 
starved since in that case the required current (j') and the charge density (n).) both decline with radius as r~^. 
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